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The improved quark mass density-dependent model (IQMDD) based on soliton bag model is 
studied at finite temperature. Appling the finite temperature field theory, the effective potential of 
the IQMDD model and the bag constant B(T) have been calculated at different temperatures. It 
is shown that there is a critical temperature Tc — HOMeV. We also calculate the soliton solutions 
of the IQMDD model at finite tmperature. It turns out that when T < Tc, there is a bag constant 
B(T) and the soliton solutions are stable. However, when T > Tc the bag constant B(T) = and 
there is no soliton solution, therefore, the confinement of quarks are removed quickly. 
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It is widely believed that the fundamental theory of strong interaction is quantum chromodynamics (QCD), which 
in principle can be applied to describe most of nuclear physics. However, due to the property of confinement and 
! asymptotic freedom, QCD can not yet be used to study low-energy nuclear physics. The challenge to nuclear physicists 
• is to find models which can bridge the gap between the fundamental theory and our wealth of knowledge about low 
\& | energy phenomenology. Some of these models have been proved to be successful in reproducing different properties 
of hadrons, nuclear matter and quark matter 0, 0, 0- Q IE El ■ The quark mass density-dependent model (QMDD) is 
one of such candidates. 

The QMDD model was first suggested by Fowler, Raha and WeinerKJ, and then employed by many authors to 
study the stability and physical properties of strange quark matter 0, UM EH • This model provides an alternative 
Qh| phenomenological description of quark confinement because in this model, the mass of quarks is given by 

Oh! B 

Q . m uA = — , (1) 

43 . n q 

B 

m s = m s0 H , (2) 

n q 



x. 

where B is the vacuum energy density and n q = n u + iid + n s o- n u , n^, n S Q represent the number density of the u, d 
and s quarks respectively. We see from Eqs.Q, @, if the distance between quarks goes to infinity, the volume of the 
system becomes infinite, n q approaches to zero and m q (q — u, d, s) approaches to infinite. This is just the confinement 
condition as that of MIT bag model 2]. However, QMDD model has its advantage since the boundary condition put 
by hand in MIT model to confine quarks is abandoned. This correspondence was confirmed by Benvenuto and Lugones 
|8| . They proved that the properties of strange matter in the QMDD model are nearly the same as those obtained in 
the MIT bag model. 

Although the QMDD model can provide a dynamical description of confinement and explain the stability and many 
other dynamical properties of strange quark matter (SQM) at zero temperature, there are many difficulties when we 
extend this model to finite temperature. They are: 1) it cannot mimic the correct "temperature T vs. density p" 
deconfinement phase diagram of QCD because the quark masses are divergent and then T tends to infinite when 
n q — > [ll|; 2) Although the basic improvement of QMDD model is the quark masses depending on density and the 
quark confinement mechanism mimics, it is still an ideal quark gas model with no quark-quark interaction included. 
However, as was shown in recent RHIC experiments, the quark-quark interaction cannot be neglected in Quark-Gluon 
plasma. In order to overcome the first difficulty, a quark mass density- and temperature dependent (QMDTD) model 
[llL H2L Ha . [LlL Ha l is proposed and it was argued that B is a function of temperature as that of Friedberg-Lee (FL) 
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model [Tiilril Il8[ . The equation of state of QMDTD model are used to study the properties of strange quark star in 
Refs.|l9 | |20 | and the results are successful. To the second difficulty, as a first step, Ref.jU has introduced a coupling 
between quark and non-linear scalar field to improve the QMDD model and have investigated the properties of nucleon 
successfully. Instead of the correspondence of MIT bag to QMDD model, a FL soliton bag is formed by quark and 
scalar field coupling for improved QMDD (IQMDD) model and the quark deconfinement phase transition can take 
place. 

This paper evolves from an attempt to extend the IQMDD model 21] to finite temperature. Since the spontaneously 
broken symmetry of the scalar field will be restored at finite temperature, the non-topological soliton tends to disappear 
when temperature T increases to Tq- If the temperature increases further, the solution becomes a damping oscillation 
and such a solution can not be taken as "soliton " solution, so the quarks are to be deconfmcd. 

By using the finite temperature field theory, we will calculate the effective potential under one-loop approximation 
and find Tc ~ HOMeV. Studying the temperature dependence of soliton solutions is our first motivation. Our second 
motivation is to find the function B(T). In QMDTD model, B(T) is an input. As an ansatz, two formulae 

Bi(T) = B [l - Ai(£-) + A 2 {^rf], 0<T<T c (3) 

with Ai,A 2 adjustable parameters, and 

B 2 (T) = B [l - 0<T<T c (4) 

were introduced in Ref.0] and Ref.^lJ respectively. However, for FL model or for the IQMDD model, B(T) can be 
calculated because the vacuum energy density B equals the difference of the value of the effective potential between 
the vacua inside and outside the soliton bag, and the vacua of the effective potential depend on temperature. 

The organization of this paper is as follows. In the following section we review the IQMDD model and its numerical 
solutions for our selected parameters. In Sec. Ill, we give detailed calculation of the effective potential of the IQMDD 
model and the bag constant as functions of temperature. The soltion solutions of the IQMDD model at different 
temperatures are presented in Sec. IV, while in the last section we present our summary and discussion. 



II. THE IQMDD MODEL 

Since the details of IQMDD model can be found in Ref.J2j|, hereafter we only write down the essential formulae 
which are necessary for our further discussion. The effective Lagrangian density of the IQMDD model is |2l| 

C = ip(iB - go - m q )ijj + -d^crd^a - U(a), (5) 

which describes the interaction of the spin-^ quark fields ip and the phenomenological scalar field a with the coupling 
constant g. m q is the mass of u(d) quark, which is given by Eq.Q. In this paper we do not discuss the strange quark. 
The potential for the a field is chosen as 

U{o) = ^° 2 + ^ + ^ + B, (6) 



b 2 > 3ac 



(7) 



The condition Q ensures that the absolute minimum of U(cr) is at a — a v =/= 0. 
U (a) has two minima: one is the absolute minimum at a large value of the a field 



3|b| 
2c 



Sac 
36 1 



(8) 



with U(a v ) = 0, another is at do = 0. The former corresponds to the physical or nonperturbative vacuum, and a = a v 
representing some condensates. The other represents a metastable vacuum where the condensates vanishes, with an 
energy density B relative to the physical vacuum. The "bag constant" B can be expressed as 



B 



2! 



3! 



4! 



(9) 
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From Eq. J5J) , the classical Euler-Lagrangian equations can be obtained as 

(ffi - ga - m q )ip = 0, (10) 

Da+^+gW = 0. (11) 
da 

In the mean-field approximation (MFA), the scalar field a is taken as time-independent classical c-number field, and 
we only consider a fixed occupation number of valence quarks (3 quarks for nucleons, and quark-antiquark pair for 
mesons). Quantum fluctuation of the bosons and effects of the quark Dirac sea are thus to be neglected. In the 
following , we will discuss the ground state solution of the system. 

In the spherical case, the a is spherically symmetric, and valence quarks are in the lowest s-wave level. Then the 
scalar field a and the Dirac equation functions can be written as 

a(r,t) = a(r), (12) 
V>(r,i) = e-"'5> 2 , (13) 



where the quark Dirac spinors have the form 



By using Eas. lfTU|) - lfH|) . we obtain 



—^— = -(e + m q +ga(r))v(r), (15) 
dv(r) 2 

—v(r) + (e — m q — ga(r))u(r), (16) 



dr r 
d 2 a(r) 2 da(r) dU 



Ng(u 2 (r) ~v 2 (r)). (17) 



dr 2 r dr da 
The quark functions should satisfy the normalisation condition 

4tt J r 2 (u 2 (r) + v 2 (r))dr = 1. (18) 

The number of quarks would be N — 3 for baryons and N — 2 for mesons. In the following discussions, we only 
constrain in the case N — 3. These equations are subject to the boundary conditions which follow from the requirement 
of finite energy: 

m- 0.^-0, 

dr 

u(oo) = 0, ct(oo) = a v . 

If we consider quarks in the lowest mode with energy e, the total energy of the system is 



E = Ne + 4tt / r 2 



1 fda x 2 



2 V dr 



U(a) 



dr. (19) 



This model has four adjustable parameter g, a, 6, c which can be chosen to fit various baryon properties such as 
the proton change radius r cp , the proton magnetic moment fi p and the ratio of axial- vector to vector coupling gA/gv 
In Ref.[2l]], a wide rage of parameters has been used to calculate above quantities. Hereafter we take one set of 
parameters a — 70.0/m -2 , b = — 2201.8/m~ 1 , c = 20000 and g — 9.8 to study the temperature dependence of the 
soliton solution. It has been proven in Ref. j2J| that this set of parameters can describe the properties of nucleon at 
zero temperature successfully. 

In Fig^we plot the a and quark fields profiles in arbitrary unit as functions of r for the above parameters. The 
quark density u 2 (r) — v 2 (r) distribution versus radius are plotted in Fig|3 
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FIG. 1: The a and quark fields in arbitrary unit as functions of r for the parameters taken as a = 70.0/m 2 , b - 
c = 20000 and g = 9.8. 



-2201.8/m -1 , 



CO 
i— 

CD 



C/3 

c 

CD 
T3 

CD 
3 
D" 




FIG. 2: The quark density u 2 (r) — v 2 (r) in arbitrary unit as functions of r for the parameters taken as a = 70.0 fm 
b = -2201.8/m -1 , c = 20000 and g = 9.8. 



III. EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 



The appropriate framework for studying phase transitions is finite temperature field theory or thermal field theory 
Within this framework the finite temperature effective potential is an important and useful theoretical tool. The idea 
of such techniques could trace back to the 1970's when Kirzhnits and Linde first proposed that symmetries broken 
at zero temperature could be restored at finite temperature)^ [23|- Subsequently Weinberg]^, Dolan and Jackiw 
25]as well as many others have adopted the effective potential as the basic tool to study the symmetry-breaking phase 
transitions. In this section we briefly outline the relevant results of Dolan and Jackiw|2fij. since the 1QMDD model is 
very similar to the model studied by those authors. 

The effective potential can be calculated to one-loop order by using the methods of Dolan and Jackiw |25j. But it 
is pointed out in Ref|16ll27| that, as an approximation, all a quantum loop diagrams may be ignored due to the fact 
that a is only a phenomenological field describing the long-range collective effects of QCD, its short-wave components 
do not exist in reality. Therefore for the rest of this discussion we shall ignore quantum corrections and concentrate 
on those induced by finite temperature effects. Then the one-loop contribution to the effective potential is of the form 



V{a-f3)=U(a) + V B (a;P) + V F (a^), 



(20) 
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FIG. 3: The temperature-dependence of the one-loop effective potential V(a;f3). At low temperature, there are two minima 
and one of which is going to disappear when the temperature reaching a high temperature T = 149.7MeV. The critical 
temperatures is set at Tc — HOMeV when two minima are euqal. 



where U(a) is the classical potential of the Lagrangian. Vp (<j\ 0) and Vp{a; 0) are the finite temperature contributions 
from boson and fermion one-loop diagrams respectivelv|25| |26j . It should be pointed out that, in this approximate 
method for the calculation of the effective potential, the plane wave quark states have been used, while, in fact, 
the quarks are confined in the finite-size solitonic configuration. Therefore, this approximate method is not fully 
self-consistent. For simplicity, we have substituted c(T) for <r(T). These contribute the following terms in the 
potential p5l| 



V b (<t;0) = 



1 



2tt 2 4 



dxx 2 \n ( 1 - e -(» a +0'O 



(21) 



i r~ 



dxx In ( 1 



-(x 2 +/3 2 m^„) 



(22) 



where the minus sign is the consequence of Fermi-Dirac statistics. m a and m q are the effective masses of the scalar 
field a and the quark field, respectively: 



B(T) 



= a + ba(T) + -a 2 (T). 



(23) 
(24) 



We fix rr? a by taking its value at the physical vacuum state [2^. In this work, B(T) is defined as the difference 
between the vacua of the effective potential inside and outside the solion bag. This means that for T < Tc, B(T) 
is the difference between the effective potential values at the perturbative vacuum state and values at the physical 
vacuum state 



B(T) = V(a Q ;0)-V(a v ;0). 



(25) 



For T > Tc, B(T) = due the fact that the vacua inside and outside the soliton bag are equal. This will be analyzed 
in more details in next section. 

From Eas. (|20|) . (|23Jl and l)24[l. we can numerically solve the effective potential V(a; 0) for different temperatures. In 
FigCH we plot the one-loop effective potential V(a; 0) as a function of a at T = lOMeV, T = 80MeV, T = HO.OMeV, 
T = 149.7McV and T = 160MeV. It can be seen from FigOthat there exist two particular temperatures. One is that 
the effective potential exhibits two degenerate minima at Tc — HOMeV which is defined as critical temperature, the 
other is that the second minimum of the potential at a ~ <r v disappears at a higher temperature T ~ 149.7MeV. 

For low temperatures the absolute minimum of V(a;0) lies close to a v , and there is another minimum at <7q. 
The physical vacuum state at a v is stable and correspondingly quarks are in confinement. As the temperature 
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FIG. 4: The bag constant B(T) as functions of T. 
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FIG. 5: The bag constant B(T), B 2 {T) and B 3 (T)as functions of T. 



increases the second minimum of the potential at Co decreases relative to the first one. At the critical temperature 
Tc — HOMeV, the potentials at the two minima are equal. The physical vacuum becomes unstable. Since B(T) = 
when T is lager than Tc, there is no bag constant to provide the dynamical mechanism to confine the quarks in the 
bag. Moreover, there is no soliton-like solution in the model when the temperature is above Tc due to the effective 
potential and vacuum structure, there only exists the damping oscillation solution and such a solution can not produce 
a mechanism to confine the quarks in a small region, we will investigate the damping oscillation in next section more 
detail. Therefore as the temperature is above Tc, the confinement of the quarks is removed completely. 

Using the obtained temperature-dependent effective potentials, from Eq.JSBJ) we illustrate the temperature depen- 
dence of the bag constant B(T) in Fig01 and it is shown that the bag constant decreases continuously with increasing 
temperature. At the critical temperature Tc, B(Tc) = 0. We can fit the analytical formula of the bag constant as a 
function of T 



B 3 (T) = B Q 



1 - cl 



T 
Tc~ 



c2 



T 



(26) 



where Bo is the bag constant at zero temperature, the parameters are cl = —0.113 and c2 = 1.113 respectively. 
For comparison, we plot B(T), B 2 (T) of Eq.QJ and B 3 (T) in Fig|5] B 2 (T) decreases smoothly with increasing 
temperature, while B(T) and B^{T) decrease dramatically with increasing temperature when the temperature is 
around 85MeV. From Eq.©, we can give the parameters A% ~ —0.67 and A 2 ~ 1.55. 
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FIG. 6: The soliton solutions for different temperatures when T < Tc- 



IV. THE SOLITON SOLUTIONS AND DECONFINEMENT PHASE TRANSTIONS 



Soliton models such as those described in Section II can be modified so as to allow for the effect of thermal 
background with temperature (3, just by replacing the relevant classical potential function U(a) by an appropriately 
modified temperature-dependent effective potential V(a;(3) 29] 30]. 

For finite temperature soliton solutions will take the same form as the one at zero temperature. But the equation 
of motion for the a field should be replaced as 

dMr)_ + 2 dap_ dV(0_ = _ 
dr z r dr da 

where the temperature-dependent effective potential V{a; (3) is defined in Ea. (|20|l . And the quark functions should 
also satisfy the normalisation condition 

4tt J r 2 (u 2 (r) + v 2 {r))dr = 1 (28) 

The functions cr(r), u(r) and v(r) also satisfy the boundary conditions following from the requirement of finite 
energy: 

«(0)= 0,^ = 0, (29) 
u(oo) = 0. (30) 

However, the situation is changed for er(r) as r — > oo. When T < Tc, in order to satisfy the requirement of finite 
energy of the solition (or other topological defects), as r — > oo, cr(r) should be equal to a v , where the potential V(a) 
has an absolute minimum. For example, at zero temperature, we take the asymptotic value (vacuum values) a — > a v 
as r — > oo, while for finite temperature, a — > <r v (j3) as r — > oo, because a v is a temperature-dependent function, and 
of course, V(cr„(/?); (3) has an absolute minimum. For T > Tc, the physical vacuum becomes unstable, and the stable 
vacuum is the perturbative vacuum which is the absolute minimum of the effective potential. Therefore, to satisfy the 
requirement of finite energy of the solition at T > Tc, we should take the asymptotic value (vacuum values) a — > 
as r —> oo. Based on above analysis, one obtain the following boundary condition for the function cr(r) as r — * oo: 

o = a v (f3), for T<T C , (31) 
a = a , for T c < T. (32) 

The set of coupled nonlinear differential equations l|15(l , l|16(l and l|27|l can be solved numerically with the boundary 
conditions Eas, 129|l . (|30|) and l|31|) when the temperature is T < Tc- In Fig[f)| we plot the soliton solutions by taking 
the temperatures as T — OMeV, 60MeV, 90MeV and HOMeV. From FigEl we can see that with the temperature 
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FIG. 7: Radius as functions of temperature when T < Tc- 
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FIG. 8: The solutions for the set of coupled nonlinear differential equations 115H . 1161 and 1271 for different temperatures when 
T c <T. 

increasing. cr v (/3) is nearly constant near a,., while <r(0) is changed dramatically. Similarly to the discussions of 
Ref.|l5l|. where they predicted that the radius of a stable strangelet increases as temperature increases, we can give 
an exact variation of radius as functions of temperature in Fig[7| which reveals that the radius of the soliton bag does 
increase with increasing temperature. 

Similarly, the set of coupled nonlinear differential equations (|15fl . (|l(jfl and l|27|l can be solved numerically with 
the boundary conditions Eas. (|29(I . H3()|l and (|32(l when the temperature is Tc < T. In FigEl we plot the solutions by 
taking the temperatures as T — 115MeV, 130MeV and 140MeV. When the temperature is higher than the critical 
temperature Tc, the bag constant is zero, the physical vacuum becomes unstable, and the perturbative vacuum state 
is stable, then the shapes of solutions are very different from that of T < Tc- With the increase of temperature, the 
soliton solutions tend to disappear. From Fig[H]we can see that, unlike the conventional soliton solutions ploted in 
FigO here the solutions have the damping oscillation, and we can not find soliton solution anymore. 

For investigating the stability of solion at different temperatures, it is necessary to calculate the total energy of the 
system at various temperatures by using Ea. (|19|) . accordingly the relevant classical potential function U(a) should 
be replaced by the temperature-dependent effective potential V{a;f3). In Fig|§l We plot the total energy of the 
system (soliton mass or energy) as functions of temperature for < T < Tc, whereto is the soltion energy at zero 
temperature. It can be seen that the total energy of system will decrease as the temperature increases. When the 
temperature is around the critical temperature Tc, the total energy of the system will dramatically decrease to 74% of 
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FIG. 9: The total energy of system (soliton mass or energy)as functions of temperature when < T < Tc, where-Eo is the 
soltion energy at zero temperature. 

Eq. Since there is no soltion solution when T > Tc, the soltion mass (or energy) for T > Tc does not exist. However, 
from the Fig|S]and Ea. ((T§|> . we can find that the total energy of system is almost equal to the quark field energy 
this also implies that the soltion bag is melted out and has disappeared when the temperature is above the critical 
temperature Tc- Since there is no soliton bag when the temperature is above the critical temperature Tc, there is no 
mean-square charge radius for the proton for T > Tc- 

As mentioned in Sect. Ill, B(T) is defined as the difference between the vacua inside and outside the solion bag. 
For T < Tc, the vacuum inside the soliton bag is the metastable vacuum, the vacuum outside the soliton bag is 
the real physical vacuum. This can be seen in Fig|S] As the temperature is above the critical temperature Tc, the 
metastable vacuum becomes the absolute one. It can be seen from FiglSI the behavior of a(r) changes dramatically 
when T > Tc- u(r) is very close to zero at any r. Therefore, the values of the effective potential are always close to 
its stable minimum at a = 0. The state at a ~ VviP) can never be realized in this case. This means that the bag 
constant B = 0. 

Based on above analysis, there is no soliton solution and the bag constant B(T) is zero for T > Tc, so there exists 
no more mechanism in soliton bag model to confine the quarks, and the quark is to be decofined when T > Te- 



la this paper, we have investigated the deconfinement phase transition of the IQMDD model at finite temperature 
and obtained the effective potential of the IQMDD model at different temperatures. We also have gotten the function 
B(T) and fitted the analytical formula of the bag constant. It is shown that the two minima of the potential at zero 
temperature will be equal at a certain temperature Tc — HOMeV. When the temperature T > Tc, the original 
perturbative vacuum state becomes stable and the original physical vacuum state becomes a metastable one. Because 
<j(r) is very close to zero at any r for T > Tc, the effective potential takes always its value at the absolute minimum 
Via ~ 0). This gives B(T) = for T > Tc- As the temperature approaches another higher temperature T ~ 150MeV, 
the effective potential has a unique minimum. Our results are qualitatively similar to that obtained in conventional 
soltion bag model at finite temperaturepMl3fll3^|. 

In IQMDD model, the confinement of the quark requires the existence of the soliton solution, and the latter depends 
on the effective potential at finite temperature. In order to investigate the behavior of the soliton solution at finite 
temperature, we numerically solve the set of coupled nonlinear differential equations. Our results show that when 
T < Tc, there exist the stable soliton solutions in IQMDD model, but when the temperature is higher than the 
critical temperature Tc — HOMeV, there is only damping oscillation solutions and no soliton-like solution exists, 
and the quarks can not be confined by such solutions anymore, then the confinement of quarks are removed and the 
deconfincd phase transition takes place. We also obtain that the radius of the soliton bag does increase when the 
temperature increases. At T — Tc the soltion bag disappears. 

Like the conventional soltion bag model, there are four adjustable parameters in the IQMDD model and the 
numerical results are also parameter-dependent. Another shortcoming of these bag models is lack of chiral symmetry. 



V. 



SUMMARY AND DISCUSSION 
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This symmetry suggests the use of a linear sigma model, which couples the quarks to pion fields as well as a scalar 
field. Solitons in such a model have been studied by Birse and Banerjee 33||3J| and Kahana, Ripka and Sonij^. It 
is of interest to extend their work to finite temperature and discuss the soliton solutions at different temperatures. 
Unlike the QMDTD model, where the bag constant B is put in by hand through an ansatz, we have obtained the bag 
constant B as a function of temperature, so our work can be used to study the effect of s quarks for strange quark 
matter. All these works are in progress. 
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